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I. INTRODUCTION 



This thesis argues that the Navy has a need for new weapons to provide ship self- 
defense against sea-skimming missiles in littoral waters. Since the Navy’s Directed 
Energy Office has already selected the free electron laser (FEL) for developmental 
funding, this thesis describes the FEL as a candidate weapon. The primary thrust of the 
thesis is to describe laser damage experiments conducted at the Thomas Jefferson 
National Accelerator Facility (TJNAF). These are the first experimental tests that study 
the damage from a short pulsed laser at a high repetition rate with a few hundred Watts of 
average power. The unique idea advanced in this thesis is that scaling rules can be 
developed that will allow these small-scale damage experiments to represent the damage 
from a large, MW-scale weapon. 

A. SHORTCOMINGS OF CURRENT SHIP SELF-DEFENSE SYSTEMS 

With the collapse of the Soviet Union and the fading from memory of the cold 
war era bipolar world, a “new world order” was supposed to have emerged. There were 
rosy predictions of a world without conflict where the nations of the earth coexisted in 
peace and prosperity. Now less than ten years after the fall of the Berlin Wall, we see that 
those predictions were grossly optimistic, and that new challenges have arisen for the art 
of diplomacy and for military science. We have had constant bloodshed from the Gulf 
War to the never-ending ethnic violence in the Balkans. In addition, many countries 
whose interests do not coincide with the interests of the United States seem to be 
vigorously pursuing programs to develop weapons of mass destruction and advanced 
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missile technology. Many of these countries’ programs seem to have been enhanced by 
the break-up of the Soviet Union, which has made weapons technology and technicians 
available to the highest bidder. 

The U.S. military by contrast has been obliged to “do more with less” due to 
budget cuts, while changing weapons procurement and tactics to meet emerging threats. 
For the U.S. Navy, there is no longer another navy in the world to pose a threat on the 
high seas, and the primary focus has shifted to the littoral with emphasis on power 
projection ashore from the sea and support of land forces. This new operating 
environment has revealed vulnerabilities of U.S. Naval forces. These vulnerabilities were 
not anticipated when current ship’s systems were being designed and built during the 
cold war’s uncertain aftermath. The primary vulnerability of U.S. Navy ships is attack by 
high-speed anti-ship missiles and is exacerbated by operating in the littoral environment. 

1. Surface-To-Air-Missiles (SAMs) 

The first line-of-defense for most ships is the Standard Missile, which is launched 
vertically and designed to destroy incoming missiles at long range. These systems, along 
with ships’ long-range radars, are ideally suited to the open ocean environment with 
visibility limited only by the curvature of the earth. In the littoral on the other hand, 
conditions are often quite different. Because ships often operate within a few miles of 
land, the detection range for incoming missiles is severely curtailed with subsequent drop 
in reaction time. Due to the speed of modem anti -ship missiles, and the short ranges of 
the littoral, the Detect-To-Engage (DTE) sequence for the Standard Missile will be 
impossible to execute in the time available. 
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Furthermore, these Standard Missile systems are not installed on all Navy ships, 
but only on the cruisers and the destroyers, about two-thirds of the ships. The rest of the 
ships have no “long range” anti-ship missile defense at all, for either the open ocean or 
the littoral. The system is very expensive both in terms of money and in terms of the 
space required. Since a ship has a very limited budget of space, and the missile magazine, 
the missile launching system, and associated radars are large and heavy, it is impractical 
to ever put them on logistics ships or amphibious ships. Cruisers and destroyers will have 
enough problems defending themselves from missile attack in the littoral, and will 
therefore be completely incapable of defending other ships. 

2. Close-In-Weapon-System (CIWS) 

All ships that might operate in the littoral do have the CIWS, which is a short- 
range gun and radar system designed to destroy an inbound missile just before it gets to 
the ship. CIWS has never been tested in battle conditions. Computer simulations indicate 
that it is unlikely the system can destroy an anti-ship missile at ranges greater than 100- 
200 meters. Although the gun has a nominal range of 2000 meters, penetrator dispersion 
seriously degrades the effective range of the system. CIWS has an average dispersion of 
two to three milliradians or 0.17 degrees, which leads to a probability of a bullet hitting 
the target as shown in Figure (1). The system has a probability of hitting the target of 
very nearly 0 % from its maximum range all the way in to about 500 meters. The 
probability does not rise above 10% until the inbound missile is within 200 meters of the 
ship. As the missile approaches the ship the probability of a hit increases and the 
accumulated hits can be calculated by applying the appropriate probabilities for each 1 00- 
meter increment of the missiles path. 
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Single Round Intercept Probability 
CIWS Height 25m, Dispersion 3mrads, Missile Height 5m 




Figure 1. Single Round Intercept Probability 



Assuming a relatively slow, non-maneuvering missile, and assuming it takes an 
average of eight hits to kill a missile, Figure (2) shows that the kill range will be 100-200 
meters as mentioned above. Killing the missile with multiple hits is not the end of the 
story. Once the missile breaks up, the fragments can still hit the ship and cause damage. 
The computer simulation was extended to analyze the trajectories of post-kill missile 
fragments. 
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Cumulative Hits 

CIWS Height 25m, Dispersion 3mrads, Missile Height 5m, 4500 rounds/min Firing Rate 




Range (m) 



Figure 2. Cumulative Hits vs. Range 

The shorter the kill range, the more likely that many missile fragments will strike the 
ship. Figure (3) shows this to be the case, and in fact at the likely kill range for CIWS, 
about 50% of the missile fragments will hit the ship with the average fragment having a 
mass of 40Kilograms, a velocity of 50 meters/second and kinetic energy of 50000 Joules. 
If the missile could be destroyed at 1000 meters or greater from the ship the probability 
of damage from missile fragments would be small. CIWS, however, is not capable of 
killing targets at those ranges. 



5 



Probability of a Fragment Striking the Ship 
Missile Speed 600m/s, Missile Height 5m, 10000 Fragments Sampled Every 100m 




Figure 3. Probability of Fragment Striking Ship vs. Kill Range 

So even under ideal conditions for CIWS with ample warning to commence 
engagement at maximum range, a slow, non-maneuvering missile leaves the phalanx 
system inadequate for ship self-defense. Under more realistic littoral conditions with little 
or no advance warning of an attack, and a fast, maneuvering missile the CIWS will be 
nearly useless. A new ship self-defense system is required. 

B. DIRECTED ENERGY WEAPONS 

One promising possibility is to use beams of high energy focused on an incoming 
missile to destroy it at long range. A major advantage of this type of weapon is that it 
travels at 300 million meters per second. Since the missile can only travel at a few 
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hundred meters per second, any maneuver the missile makes becomes meaningless. 

There are three major methods of generating the large energies necessary to destroy an 
anti-ship missile in flight: Particle Beams, High Powered Microwaves (HPM), and 
Lasers. All of these technologies have been studied for many decades to determine their 
suitability as anti-missile weapons. The advantages and disadvantages of each will be 
discussed in turn to see where the scientific research has lead so far. 

1. Particle Beams 

The idea of using a high-energy beam of subatomic particles, electrons, protons or 
neutrons, to shoot down a missile has been around since the early 1970’s, but was given 
major consideration and funding as part of the Strategic Defense Initiative (SDI) in the 
1980’s. Of all the possible Directed Energy Weapons the Particle Beam seems to face the 
most serious obstacles. 

The main problem for an electron beam anti-missile weapon is propagation loss, 
which has three major components: Ionization Loss, Bremstrahlung Loss, and Beam 
Divergence due to elastic scattering. Ionization Loss occurs when energy from the beam 
ionizes surrounding air molecules and amounts to approximately 0.34 MeV/m. This 
means that a 500 MeV electron beam would travel less than 1500 meters even if there 
were no other loss mechanism. But there are other loss mechanisms. Bremstrahlung Loss 
is caused when beam electrons are accelerated in curved paths around atomic nuclei due 
to coulomb interaction. The electron beam then loses energy by emitting Bremstrahlung 
photons. Bremstrahlung Loss alone would limit a 500 MeV beam to a propagation path 
of about 360 meters. Elastic scattering of the beam occurs when electrons “collide” with 
more massive air molecules as in Bremstrahlung Loss and the trajectories of the electrons 
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are changed. This causes the electron beam to spread as it propagates through the air. 
Elastic scattering alone would cause a 500 MeV electron beam to dissipate in 240 meters 
[Ref. 1]. Actual experimentation with all loss mechanisms working together shows that 
an electron beam only propagates a few meters through the air. 

2. High Powered Microwaves (HPM) 

HPM technology seems to hold more promise than particle beams because of the 
maturity of the field and we know we can propagate them through the air. But HPM 
weapons also face challenges. Assuming that an adequate source of HPM is available 
there are generally two methods of employment to defeat incoming missiles. One is to 
aim an intense beam at a specific target to destroy it at long range, and the other is to 
sweep a large area in the hope of disabling many targets. The second method is 
unattractive for military applications because of the lower probability of kill for each 
target and because of the increased likelihood of fratricide. 

Once a source is found and a tactic is chosen, the HPM radiation must reach the 
target by propagation through the atmosphere. Atmospheric attenuation of HPM at sea is 
dominated by suspended particles such as dust and water droplets, which condense 
during rain or fog. Figure (4) shows that in the absence of fog or rain, a water vapor 
absorption peak at 22 GHz and an oxygen absorption peak at 60 GHz dominate 
attenuation of HPM. There is a relative minimum in attenuation at 35 GHz. Once the 
radiation reaches the target there are two basic transmission paths into the electronics of 
the missile. The “front-door” path is one designed to transmit microwaves in the normal 
operation of the system, and the “back-door” path is one of radiation coupling through a 
path not designed for transmission. It requires less energy on target to defeat a missile by 
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the front-door path, but there is certainly no guarantee that the enemy will cooperate and 
design missiles appropriately. 



ATTENUATION P-WAV, <«/»•) 




FREQUENCY (CNtJ 



Figure 4. Atmospheric Attenuation of HPM, From Ref. [2] 



An additional problem for HPM weapons is the ease of shielding against them. 
Herbert, using a digital watch and a microwave oven, conducted simple experiments that 
were quite enlightening. A three-second exposure to microwaves caused an unshielded 
watch to quit working. When a similar watch was placed in a box made of aluminum foil 
and placed in the microwave, it was able to withstand repeated exposure including 
individual exposures of as long as 90 seconds with no loss of function [Ref 3]. Assuming 
a relatively slow missile capable of 300 m/s, the shielded missile would travel 27 
Kilometers during that 90 seconds and the ship would be dead. 
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Another problem for long wavelength HPM is diffraction. Beam spread angle is 
proportional to wavelength by the relation 0 = 1 .22 A, / D, where D is the aperture size. 
Beam spot area is then A = n (R0) 2 , where R is the range to the target. With an antenna 
diameter of 1 0 meters and a wavelength of 6 cm (5 GHz), the beam spot is 4208 m 2 at a 
range of 5 km. The energy is spread over a wide area and the effectiveness of the weapon 
reduced. 

3. Lasers 

Lasers are in use in our society today for so many mundane purposes that we take 
them for granted. We use them as pointers for presentations in offices, to play our music 
and movies at home, to cut things in industry, and to perform surgery in hospitals. Laser 
weapons have excited the American imagination for decades. In the 1950's early science 
fiction shows had "ray guns", in the 1960's Star Trek had "phasors", in the 1970's and 
1980's Star Wars had "Death Stars" with lasers which could destroy entire planets at a 
shot. The scientific community has worked hard in conjunction with the military to 
explore the potential of the laser as a weapon. The first proof of this potential was the Air 
Force's ALL project in the 1970's, followed by the Navy's MIRACL project in the 1980's, 
and the Air Force's ABL program in the 1990's [Ref. 4], There are many different types 
of lasers that have different strengths and weaknesses, and all must be judged against the 
stringent requirements of a shipboard laser weapon. 

The laser must be able to emit enormous power, but be small enough to fit on a 
ship. It must operate at a wavelength that propagates well through the atmosphere, so that 
range does not suffer. It must not produce dangerous byproducts that cannot be disposed 
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of at sea. We shall look more closely at each one of these criteria and see how the 
available lasers stack up. 

4. Atmospheric Propagation 

It is important to study the transmission of laser radiation through the atmosphere, 
because certain wavelengths, and therefore certain lasers, can be eliminated on the basis 
that they cannot propagate 5 km to the target with sufficient energy to destroy a target. A 
non-linear phenomenon known as thermal blooming can occur with high power lasers 
when the air is stagnant and the target is approaching on a constant bearing. The air 
between the ship and the target is heated, changing the index of refraction, and forming 
an effective lens that can disperse the laser beam. Since thermal blooming is a non-linear 
effect, the more power in the beam, the faster thermal blooming will occur. As a result it 
is critical to find a wavelength with the minimum absorption and, therefore, maximum 
transmission through the atmosphere. Figure (5) shows atmospheric transmittance for 
wavelengths of interest for lasers, but the finer detail is necessary to pick an exact 
window for optimal transmission. There are regions of high transmittance in the 3-5 
micron band and the 8-12 micron band, which have been the typical wavelengths for IR 
sensors. There are much smaller sections of the spectrum, which give close to 100% 
transmittance at 1.042 pm, 1.06 pm, 1.6 pm, 2.2 pm, and 3.8pm. Cook and Albertine 
have done a more detailed study of this issue in relation to a maritime deployed high- 
energy laser weapon system (HELWS). Their conclusions are summarized in Figures (6) 
and (7), which show the most promising wavelength to minimize absorption and the 
thermal blooming problem is 1 .042 pm. 



ll 




Wavelength tf/m) 



Figure 5. Atmospheric Transmittance for 0-15 pm, From Ref. [5] 




Figure 6. Total Absorption 1/km, From Ref. [6] 
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Figure 7. Total Extinction 1/km From Ref. [6] 



The FEL, which can be designed for a wide range of wavelengths, is the only 
laser capable of producing this exact wavelength. In addition to design flexibility, once 
the FEL is built it can be easily modified by changing the electron beam energy to 
operate over a range of wavelengths within a factor of two. Other lasers such as chemical 
lasers, gas discharge lasers, excimer lasers, and x-ray lasers are confined to a specific 
wavelength or wavelengths by their generation mechanism. So, an analysis of the first 
requirement, atmospheric propagation, shows the FEL to be a viable option for a laser 
weapon to defend against anti-ship missiles. 

5. Size and Power 

Experiments in conjunction with the MIRACL program indicate that an energy 
flux or power density, O ~ 10 kW/cm 2 is needed to destroy a missile with a dwell time of 
a few seconds. 
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( 1 ) 



A n(RQ) 2 ’ 

where P is power, A is spot size, R is the range, and 0 is the beam half angle in the far 
field. So for a spot size of 100 cm , a Megawatt-class laser is required. By integrating O 
with respect to time we can find the energy deposited at the target in a given time, or 
fluence, F. 

F = jW = Or , (2) 

where t is the pulse length and <I> is assumed constant over the pulse. By substituting Eq. 
(1) for <t> and the diffraction limited 0 = 1.22 A./D into Eq. (2) we get an expression for 
the fluence. 



F = 0.21 



PD\ 

m 2 * 



( 3 ) 



where X is the wavelength and D is the diameter of the aperture. Because of atmospheric 
extinction more energy is necessary at the ship to obtain the required energy at the target 
[Ref. 7]. 



F - nm ‘W T ’ (4) 

where T = e - ^ is the atmospheric transmittance, and a is the extinction coefficient, F 
[kJ/cm 2 ], P [MW], D [m], x [sec], A.[|am], and R [km]. 
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II. FEL THEORY 



A. BASIC OPERATION 

A free electron laser (FEL) takes advantage of a simple physical characteristic. 
When bunched electrons are accelerated they emit photons coherently. A static magnetic 
field, transverse to the direction of a beam of relativistic electrons, causes the electrons to 
change directions or “wiggle,” from side to side. As the bunched electrons wiggle they 
give off coherent radiation with proper design of the electron beam, the optical cavity, 
and the magnetic field. The major components of an FEL are a source of relativistic 
electrons, an undulator and appropriate optics for a resonant cavity [Ref. 8]. 

1. Electron Beam 

Electrons are produced by an electron gun, then accelerated to relativistic 
energies. The electron energy is, E = ymc , where y is the Lorentz factor, m is the 
electron mass, and c is the speed of light. Electron energies may be from a few MeV to a 
few GeV depending on the desired laser wavelength. Bunching of the electrons is 
required for optical gain to occur. The size of the accelerator is a concern for shipboard 
application of the FEL and electron beam quality is also important. Electron beam 
quality is a measure of the spread in energies of electrons throughout the beam, and is 
best described (smaller is higher quality) by the ratio of the spread of the Lorentz factors 
to the average Lorentz factor, Ay / y. A typical layout for an FEL with a recirculating 
electron beam is shown in Figure (8). 
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Figure 8. FEL Electron Source 



2. Undulator 

In the undulator, the optical field is formed and amplified by the interaction of the 
relativistic electron beam with the spatially periodic magnetic field. The magnetic field is 
produced by the alternating polarity arrangement of the magnets. A picture of an 
undulator, sometimes referred to as a “wiggler,” is provided in Figure (9). The undulator 
wavelength, Xo, is the distance along the beam axis between magnet pairs, and is given 
by, Xo = L/N, where L is the total length of the undulator and N is the number of 
undulator periods. The magnetic field can be produced by electromagnets, but most FEL 
undulators are generally constructed with permanent magnets as in Figure (9). 




Figure 9. FEL Undulator 
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3. Optical Cavity 

The optical cavity is oriented along the undulator axis and extends beyond the 



undulator length. In Figure (10), the arrows in the center represent the undulator periods, 
which are bracketed by mirrors containing the optical field. During successive passes 
through the undulator, the optical field is amplified and some fraction of the coherent 
radiation is allowed to escape and used in the weapon system. The physics of the FEL 
can be described by two equations - one for the electron motion and one for the optical 
field evolution. These equations will be described in detail in the following section. 



B. PENDULUM EQUATION 

It will be shown that the pendulum equation describes the transfer of momentum 
and energy between a free electron and the electromagnetic wave in the magnetic field of 
the undulator. The governing force equations for the electron are the relativistic Lorentz 
equations, 




FREE ELECTRON LASER RESEARCH 



Figure 10. FEL Optical Cavity 




(5) 




( 6 ) 



y- 2 =l-p.|3, 



(7) 
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where |3 is the velocity of the electron as a fraction of the speed of light. The electrical 
field and magnetic field of the optical wave are given by, 

E r = E 0 [ cosCnsinOnO], (8) 

5 r = 5 0 [sin(ncos(^),0], (9) 

where T = kz - at + <j> , the wave number is k = 2n/X, co is the angular frequency, <j> is 
the optical phase angle, and X is the optical wavelength. The undulator static magnetic 
field is represented by, 

B= fl 0 [cos(V)»sin(* 0 r),0], ( 10 ) 

where the undulator wave number is k 0 = 2n/\ , and Xq is the undulator wavelength. 
Substituting the fields in Eq. (8), (9), and (10) into Eq. (5), and resolving the transverse 
components gives, 

[E 0 ( 1 - P- )(cos( v F),- sin( v f / ),0) + B 0 P (- sin(£ 0 z ), cos(^ 0 z ),0)] , (1 1) 

at me 

where P x is the transverse electron velocity and p. is the electron velocity in the direction 
of the undulator axis. For a relativistic electron E 0 (\ - p.) « £ 0 /2y 2 , which can be 
neglected for y » 1 , leaving, 

= — ~ Bo P; [~ sin( A: 0 z ), cos(A: 0 z),0] . (12) 

ot me 

Eq. (12) can be integrated in time with the constant of integration set to zero indicating 
ideal injection of electrons. The result is, 

P x =- — [cos(V)> sin (£o z )>°]> ( 13 ) 

y 
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where we have defined an undulator parameter K = eB G \/2nmc 2 . The undulator 
parameter allows adjustment of the transverse electron velocity by adjusting the 
undulator magnetic field intensity or the spacing between magnets. 

Insert Eq. (8) and Eq. (13) into Eq. (6) to get, 

* dy eE 0 K 

y= — 2"COs(; +4>), (14) 

at ymc 



y eE K 
y = ~y T mc 2 C ° S ^ 

where the electron phase is C = (k + k 0 )z —<£>t . The initial phase at time t - 0 is 
Co = ( k + k 0 ) z o > but since k » k 0 we can approximate^ « kz 0 = 2nz 0 /X , which gives a 

relation between the electron phase and the optical wavelength. Insert Eq. (13) for P ± 
into Eq. (7) to get 

r" 2 =i-PP =i-Pi-fc 2 =7^-- 06) 

Differentiate both sides with respect to time to get, 



Y _ T% P. 
y 1 + K 2 ’ 



(17) 



then using Eq.(15) substitute P r =(^/(k + k 0 )c , and solve for t, , 



k l + K * (k + k 0 )ceE 0 K 
C = 2 o 2 C0S (C + < t>) 

y p : y me 



(18) 



Applying approximations based on k»k 0 , P z « 1 when y»l, and defining 



dimensionless time x = ct/L, we have the pendulum equation, 



£=|a|cos(£+4>), (19) 
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a = 



4 nNeKLE 0 
y 2 mc 2 



( 20 ) 



where \a\ is the dimensionless optical field, N is the number of undulator periods, L is the 



total length of the undulator, and (_) = d{_)/dx. 



C. OPTICAL WAVE EQUATION 

The second part of the mathematical description of an FEL is the wave equation. 
Beginning with Maxwell’s wave equation for the vector potential A such that 





1 8A 


(21) 


K = - 


"cUt’ 








B r =VxA, 


(22) 


-2 1 d 2 


4n - 




V 2 , -r 

c 2 dt 2 


A = — J ± , 

c 


(23) 



where J ± is the transverse current density and 



A(z,t ) = — E(z, t)[cos( v P),- sin(¥),0] , 
(0 



(24) 



where the optical phase is ¥ = kz - cot + <[> (z, t). Taking two spatial derivatives and two 
time derivatives of Eq. (24) yields. 



d^A 1 dE 
dz 2 k dz 



kJr Tz 



[cos(T),-sin( v P),0] 



+ i^[sin(T),cos(T),0] + if 



k dz 2 



dz 



[cos( v F),-sin('P),0] 



Ed 2 cb r -i E 

+ J [cos(T),- sin( v P),0] + - 






[- sin( v P),-cos(T / ),o], (25) 
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and 



co cTA dE 
c dt 2 dt 



'd$ 

dt 



-co 



8 2 E> 



[cos(T) -sin(T), 0 ] + -^ r [sin(T),cos(T) > 0 ] 



+ ■ 



dE 

dt 



dt 



-co 



[cos(T) ,-sin(T), 0 ] + £^[cos(T) -sin(T), 0 ] 



[-sin^X-cosC^XO]. (26) 

Next, assume the optical phases and amplitudes are slowly varying in time and 
space such that dE/dz « kE, d§/dz « £<J> , dE/dt « co E, d$/dt « co<() , and 
co = kc. Eliminate second-order terms and insert Eq. (25) and Eq. (26) into Eq. (23) to 
get 



+ E 



d<|> 

~dt 



-co 



"-2 i a 2 ' 


■ t 


'dE 1 dE' 


V 2 -— r — 7 


A * 2 


— + 


C 2 ar 2 _ 




_dz c dt _ 




+ 2E 








dz c dt _ 



[cos^), sin^XO] 



471 - 



(27) 



The current density for a single electron is J x = -ec(3 x . Substituting Eq. (13) for and 
introducing new time and space coordinates, t = ct / L, and z = z-ct, we have, 



]_dE 
L di 



[cos('f'), sin(T),0] + 2E [- sin( v P),- cos( v f / ),0] 



4neK r i 

— [cos(* 0 z),sin(* 0 r)]. (28) 

The dimensionless time x progresses from zero to one from the beginning of the 
undulator to the end, and the coordinate z follows a point on the optical pulse traveling 
at speed c. Insert Eq. (28) into Eq. (27), take the average (...) over many sample 
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electrons, and multiply by the electron particle density, p, giving the equations of motion 
for the optical field amplitude and phase, 



dE 

dx 

3<J) 

Ik 



2nep LK 

y 



(cos(^ +<)>)>, 



InepLK (sin(^ +<())), 



(29) 

(30) 



where C, = (k + k 0 )z - cons the electron phase. Use the complex electric field E = E 0 e '* , 
to get 



^ = -^^^[(cos(<; +<(>))- z(sin(^ +<t)))]e'* - - 2%epLK ( e < ). ( 31 ) 

or y y 

Substituting the definition of the dimensionless optical field a = \a\e‘ * from Eq. (13) we 
get the final form of the wave equation 



where the dimensionless current is j = 871 2 e 2 pNK 2 L 2 /y 3 me 2 . The growth of the optical 
field, and therefore the development of laser energy, is dependent on the dimensionless 
current and the average electron phase. So we need both a current j and electron 

bunching (e ~' Cj ) ^ 0 to produce laser energy. 
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III. FEL DESIGN AND PARAMETERS 



A. TJNAF FEL 

The most powerful FEL ever operated is at the Thomas Jefferson National 
Accelerator Facility (TJNAF). In April 1999, the laser operated continuously at 710 
Watts of average power for six hours. The FEL is designed to achieve 1 kW and will 
probably soon do so. Near term modifications now in the planning stage will boost the 
power to 20 kW. Table (1) shows the parameters of the TJNAF FEL and compares them 
to the requirements for a shipboard anti-missile defense weapon [Ref. 9]. 



Parameter 


TJNAF FEL 


Weapon FEL 


Average Power 


P= 1 KW 


P= 1 MW 


Average Current 


7=5 mA 


7 = 900 mA 


Electron Energy 


ymc 2 = 48 MeV 


ymc 2 = 100 MeV 


Lorentz Factor 


y = 94 


y = 196 


Undulator Field 


B 0 = .38 T 


B 0 = .27 T 


Undulator Wavelength 


X 0 =2.7 cm 


A^= 4 cm 


Undulator Periods 


o 

II 


<N 

II 

* 


Undulator Length 


Z, = 108 cm 


Z, = 100 cm 


Undulator Parameter 


A: =0.96 


K= 1 


Rayleigh Length 


z 0 = 40 cm 


z 0 = 10 cm 


Electron Charge/Bunch 


I e /c= 60 pC 


IJc= 1800 pC 


Peak Current 


7=60 A 


7=600 A 


Electron Beam Radius 


r b = 100 pm 


r b = 300 pm 


Pulse Length 


x = 0.4 ps 


x = 3 ps 


Pulse Repetition Rate 


PRR = 18.7/37.4/74.85 MHz 


PRR = 500 MHz 


Optical Resonator Length 


L opt = 8-° 1 m 


L op ,= 10 m 


Output Coupling 


10% 


10% 


Resonator Cavity Losses 


Q = <0.5% / pass 


Q- <0.5% /pass 


Optical Wavelength 


X - 3-6 pm 


X = 1 pm 



Table 1. Comparison of TJNAF FEL with a weapon FEL 
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The big differences are increases in the (i) peak current by a factor of 10, (ii) the 
repetition rate by a factor of 7, (iii) the electron beam energy by a factor of 2, and (iv) the 
pulse length by a factor of 7. 



B. MODELING A MW FEL 

As we can see from Table (1) there is still a long way to go from the current FEL 
to a weapons application FEL. Since we do not have a MW-class FEL to perform 
experiments, we must turn to scaling to determine the effectiveness of a FEL weapon 
without spending hundreds of millions of dollars. As mentioned in Chapter I, power 
density <t> « 10 kW/cm over a spot size A « 100 cm is required to destroy a missile with 
a dwell time of a few seconds. By scaling we can get an appropriate spot size for any 
laser power. With a power density of O = 10 kW/cm 2 , a 100 W laser must use a spot size 
of 1 mm 2 , and a 1 kW laser must use a spot size of 10 mm 2 . 

Scaling will only work, however, if the heat diffusion is independent of spot size. 
Schriempf gives a detailed treatment of diffusion, but a summary is useful. The heat flow 
equation is 

pC^ = KV 2 T+^ (33) 

where p is the material density, C is the specific heat, T temperature, t time, K thermal 
conductivity and A is the heat produced per unit volume per unit time [Ref. 10]. For a 
semi-infinite solid with no phase change and material properties that are temperature- 
independent, Eq. (33) reduces to 



d 2 T }_dT___± 
dz 2 k dt K ’ 



(34) 



24 



where the thermal diffusivity is k = K/ pC. The heat produced per unit volume per unit 
time is given by 

^(z,0 = (l-9?)O(/)ae" a -', (35) 

where 9? is the reflectivity, O(t) is the power density at depth z, and a is the absorption 
coefficient. At optical wavelengths, for most solid materials, the absorption coefficient is 
very large a « 1 0 6 /cm. The exponential term approaches zero, and A is small compared 
to the time and space derivative terms. As a result we can assume A « 0 and Eq. (35) 
becomes 



The solution to Eq. (36) is 



d 2 T 1 dT 
dz 2 k dt 



T(z,t ) 



O 0 _D 
K 




(36) 



(37) 



where <E>o is the constant power density on the target and the characteristic thermal 

diffusion length is D = 2vk7 which represents the distance required for T to drop to Me 
times its initial value. In the semi-infinite approximation we ignored radial heat flow, but 
in order for this to be valid the spot size must be much larger than D, or the target 
diameter d < D. If these conditions are not met, we expect heat to diffuse outside of the 
laser spot and the spot will not be heated effectively. 

One of the targets to be irradiated was Al-6061 . A calculation of the thermal 
diffusion length associated with heating the sample to melting temperature was 
performed with p = 2700 Kg/m 3 , C = 896 J/Kg-K, K = 180 W/m-K, k = 7.44 x 10’ 5 m 2 /s, 
T m = 855K, 0o = 10 8 W/m 2 and 
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D = 2^/i<x", 



(38) 



7tK 2 Ar 2 

40 2 k 



(39) 



Where x is the time required to bring the material from ambient temperature to melting 
temperature and AT is the temperature change. The result for Al-6061 was D = 1.79 mm. 
In order to melt through an aluminum sample, the scaled laser spot must have an area 
A » 7i(l .79 mm) 2 = 10 mm 2 or the target itself made small with d * D . These 
calculations were experimentally verified with samples of Al-6061. Using d = 1 cm, the 
sample was irradiated with no melting after several minutes. Other samples with 
d = 2 mm were melted in a few seconds. 

Another target to be irradiated was Slip-cast Fused Silica (Si02). A calculation of 
the thermal diffusion length associated with heating the sample to its melting temperature 
was performed with p = 2200 Kg/m 3 , C = 920 J/Kg-K, K = 1.26 W/m-K, k = 5 x 10' 7 
m 2 /s, T m = 1980K [Ref. 11], (To = 1 0 8 W/m 2 and the result was D = 0.021 mm. Therefore, 
with the insulating material, fused silica, the scaled laser spot must have an area 
A » 7i(0.021 mm) 2 = 0.0014 mm 2 , which is attained with the 1 mm 2 beam used at 
TJNAF. 



C. PULSE TRAIN 

Since there is growing interest in lasers for many scientific and engineering 
applications, research into short pulse effects has increased sharply in the last few years. 
The pulse train of an FEL is different from any other laser and its interaction with matter 
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at high-average power has not been studied. The TJNAF FEL has a pulse length of 
x = 0.4 ps and a repetition period of T = 27 ns illustrated in Figure (11). The duty cycle 
D is the fraction of time the laser is actually irradiating the target, 



x 4 x 10~ 13 s 

D = - = = 1-5 x 10~ s . 



T 2.7 x 10 5 



(40) 



The peak power in each micropulse P is 



- P 7\0W 

f =FIrf= 4M ' 



(41) 



Where P is the average power of 710 W. 




Figure 11. FEL Pulse Format 

Comparing the TJNAF FEL to another short pulse laser is instructive. The 
Lawrence Livermore National Labs (LLNL) 1.053-pm Ti:sapphire CPA system [Ref. 12] 
has a pulse length x = 0.4 ps, but a pulse repetition rate of only 10 Hz, so the period is 
T = 0. 1 s, and the peak power is P = 2.5TW . The duty factor is 



^ 4 x 10” 3 s „ 

D = — — = 4 x 10“ 2 , 

0.1s 



( 42 ) 



so that the average power is 
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P =(2.5x10 I2 f)(4x10‘ 12 ) = 10^. (43) 

Note that the LLNL laser has a much higher peak power than the TJNAF FEL, but the 
TJNAF FEL has more than seventy times the average power because of its high duty 



cycle. The experiments detailed in the following section were conducted to study the 
effects of the unique FEL pulse format in laser-matter interaction on small samples. 
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IV. FEL EXPERIMENTS 



A. LABORATORY LAYOUT 

Experiments for this thesis were conducted at the Thomas Jefferson National 
Accelerator Facility (TJNAF), which is operated by the Southeastern Universities 
Research Association (SURA) for the U. S. Department of Energy (DOE). Once a laser 
beam has been created, as described in chapter two, the beam is sent by a low-loss optical 
path to user laboratories in another part of the building. User laboratory number one was 
used for all experiments described below. M. Shinn, S. Benson, B. Yunn, G. Neil, K. 
Jordan and J. Gubeli operated the laser and other laboratory equipment. The optical 
bench was configured with a focusing lens, a sample holder and two video cameras. One 
camera was focused on the front of the sample and one on the back of the sample as 
shown in Figures (12) and (13). 




Figure 12. Back view of sample in user lab arrangement 
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Figure 13 Front view of sample in user lab arrangement 



B. EXPERIMENTAL PROCEDURE 

Samples were irradiated by a laser beam with wavelength X = 4.825 pm through a 
calcium fluoride lens with a focal length of 300 mm. The pulse repetition frequency 
(PRF) was 37.425 MHz and the average power meter in the user lab read 100-103 W 
with an error of +/- 5 W. Since a lens focused the beam, the beam area decreased with 
distance along the direction of propagation to a minimum waist radius of wo = 80 pm at 
the focal point. Paraxia, a beam propagation code, was used to model the beam 
diffraction and find the target position giving the desired intensity of 10 kW/cm [Ref. 
13]. Figure (14) shows a graph of irradiance versus distance from the focal point, with 
the negative numbers indicating positions in front of the focus. 
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Distance from focus (mm) 

Figure 14. Irradiance vs. Distance from Focal Point, From Ref. [13] 

An inadiance of 10 kW/cm 2 occurs when the sample position is between 25-30 mm in 
front of the focal point. The samples were actually placed 26 mm in front of the focus. 
The bum-through time was determined by observing a signal on a power meter placed 
behind the samples, and by watching for the presence of coherent harmonics in the 
visible spectrum on an iris placed approximately 15 cm behind the samples. 

C. DESCRIPTION OF RESULTS 

The Naval Research Laboratory (NRL) provided the sample materials, which 
included Slip-cast Fused Silica (Si02), Polyimide Fiberglass, and F2 Epoxy. The last two 
were from the same batch as the material used for previous continuous wave (CW) laser 
penetration studies. 
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1. Sample #1 Slip-cast Fused Silica 

The sample provided by NRL is 6.9 cm by 7.4 cm and varies in thickness from 
0.9 cm to 1 .9 cm. The front of the sample after irradiation is shown in Figure (15). 




Figure 15. Slip-cast fused silica sample 

The sample was irradiated eight times and the results are summarized in Table (2). 



Run 


Average 


Exposure 


Penetration Depth 


Number 


Intensity (kW/cm 2 ) 


Time (s) 


Rate (mm/s) 


1 


9 


9 


0.20 


2 


9 


110 


0.081 


3 


9 


13 


0.20 


4 


9 


24 


0.125 


5 


9 


41 


0.081 


6 


9 


2 


0.35 


7 


500 


3 


3.0 


8 


500 


11 


3.0 



Table 2. Irradiation of Slip-cast Fused Silica, After Ref. [14] 
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The irradiations were done from left to right. The last two runs were done with 
the sample at the focus instead of 26 mm in front so that the beam waste radius was w = 
wo = 80 pm and the intensity was 500 kW/cm 2 . The last two runs were conducted to 
investigate the effects of much higher power density. As shown in Figure (15), the first 
six runs were along the top of the sample and the last two were approximately 1 cm 
below. The second, seventh and eighth runs penetrated the entire 0.9-cm thickness of the 
Fused Silica material. Figure (16) shows the effects of exposure time on penetration 
depth rate. 



Penetration Depth Rate vs Exposure Time 




Figure 16. Exposure Time vs. Penetration Depth Rate for Fused Silica 
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As the exposure time was increased more smoke and debris filled the hole blocking the 
path of the laser beam and causing the penetration depth rate to decline over time. In the 
future, we can explore whether altering the FEL wavelength during sample irradiation 
improves penetration depth rates through smoke and debris. 

Since the second run was the only one to punch completely through the material 
that was also conducted at the primary power density of interest, we shall take a closer 
look at it. A digital picture of run two was taken through an optical microscope as shown 
in Figure (17). 




Figure 17. Close-up of damage to slip-cast fused silica in run 2 

Although the beam diameter was only 1 . 1 mm, the melted portion at the surface of the 
sample measured 5 mm in diameter. The hole is tapered with the melted portion on the 
back of the sample measuring only 2 mm in diameter. The reasons for this are twofold. 
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First, the beam profile is Gaussian in nature with the highest intensity in the center of the 
beam and intensity down by 1/e at the beam radius. Second, the sample face is 26 mm in 
front of the focus and the sample back is 1 7mm from the focus, so the beam size is 
decreasing as it proceeds through the material as shown in Figure (18). 



1.2 





Figure 18. Sketch of beam focusing effect 

The vertical scale in Figure (18) is exaggerated by a factor of five with respect to the 
horizontal scale in order to demonstrate the effect. 

Examination of the hole from run two through an optical microscope reveals a 1- 
mm thick layer of melted, and rehardened, SiC >2 filling the hole at the back of the sample. 
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It was clear from the video and the rear power meter that bum-through occurred in run 
two, but melted material solidified and sealed the hole at the back of the sample. A 
picture of the back of the target taken through a Scanning Electron Microscope (SEM), 
Figure (19), shows the hole from run seven is fairly irregular with a great deal of debris. 




Figure 19. SEM photograph of damage to Fused Silica sample from run 7 

The volume of the hole in run two is estimated by, 

V= n[R(z)f dz , (44) 

where the radius approximately changes linearly as 

R(z) = 0.53mm -0.02z, (45) 

which gives a volume of V— 5.6 mm 3 . By a similar calculation, the volume of the entire 
damaged region, including the melted and rehardened portion, is estimated to be 
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V= 92 mm 3 . Based on the density of fused silica of p = 2.2 g/cm 2 , the amount of 
material removed was 0.012 g, and the amount of material damaged was 0.20 g. The heat 
energy deposited during run two is given by, 

E = Px = OAx , (46) 

which gives E = 9.7 kJ deposited during the 110 second run. The heat of ablation is then 
48 kJ/g for the damaged area. In future studies this value should be correlated to the 
values obtained for similar materials in much larger scale CW laser experiments in the 
past, with appropriate compensation for experimental environment and wavelength. 

2. Sample #2 Polyimide Fiberglass 

The sample provided by NRL is 1 1.4 cm by 10.1 cm and 2 mm thick. The 
damaged area of the sample, after irradiation, is shown in Figure (20). 




Figure 20. Polyimide Fiberglass target 

The sample was irradiated three times and the results are summarized in Table (3). 
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Run 


Average 


Exposure 


Penetration Depth 


Number 


Intensity (kW/cm 2 ) 


Time (s) 


Rate (mm/s) 


1 


9 


7 


0.28 


2 


9 


2 


0.90 


3 


9 


1 


1.1 



Table 3. Irradiation of Polyimide Fiberglass, After Ref. [14] 



The irradiations were done from left to right with the sample 26 mm in front of 
the focus of the beam. Only the first run achieved bum-through of the material, with the 
entry hole 3 mm in diameter and the exit hole 1.5 mm. All three holes show significant 
charring adds an additional term to the heat transport equation and impedes ablation. 
Investigation with an optical microscope reveals a raised lip of material around the face 
of the hole that does not appear on the fused silica sample and much more roughness as 
observed in Figure (21 )-(23). The charred region extends to a diameter of 8.0 mm for run 
one, 6.5 mm for run two, and 5.4 mm for run three. The lip height is 0.3 mm for run one, 
0.1 mm for run two, and 0.05 mm for run three. These measurements indicate that as the 
dwell time increases, the radial extent of the damage area increases, and more material is 
deposited around the edge of the hole. There is no evidence of melted and rehardened 
material present in the holes as found with the fused silica indicating a different 
mechanism for ablation in the two samples. 
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Figure 21. Close-up of damage to Polyimide fiberglass in Run 1 




Figure 22. Close-up of damage to Polyimide fiberglass in Run 2 
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Figure 23. Ciose-up of damage to Polyimide Fiberglass in Run 3 



3. Sample #3 F2 Epoxy 

The sample provided by NRL is 10.0 cm by 1 1.5 cm and 1.5-mm thick including 
a 1 ,6-cm thick polyurethane foam backing. The damaged area of the sample, after 
irradiation, is shown in Figure (22). 




Figure 24. F2 Epoxy sample 



The sample was irradiated three times and the results are summarized in Table (4). 
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Run 


Average 


Exposure 


Penetration Depth 


Number 


Intensity (kW/cm 2 ) 


Time (s) 


Rate (mm/s) 


1 


9 


11 


0.10 


2 


9 


6 


0.12 


3 


9 


3 


0.10 



Table 4. Irradiation of F2 Epoxy, After Ref. [14] 

In each case, it appears that the F2 Epoxy was completely penetrated and the 
ablation of the foam backing had begun, but not completed. The videotape showed 
flames engulfing the upper portion of the sample and Figure (24) shows the black charred 
area extending to the edge of the sample. Significant charring was evident when the 
sample was viewed with the optical microscope, very similar to the Polyimide sample. 
There was also evidence of some melting, but not as much as occurred in the Fused Silica 
sample. The holes appear to be filled with the charred debris of the polyurethane 
backing, making hole depth measurements difficult and rendering penetration depth rates 
unreliable. Figures (25)-(27) show the details of runs one through three. 




Figure 25. Close-up of damage to F2 Epoxy in run 1 
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Figure 26. Close-up of damage to F2 Epoxy in run 2 




Figure 27. Close-up of damage to F2 Epoxy in run 3 

The damaged region extends to a diameter of 1 1.3 mm for run one, 7.5 mm for run two, 
and 5.2 mm for run three. There is a lip around each of the holes, but much smaller than 
the polyimide sample showed. The lip for rim one was 0.05 mm and for runs two and 



42 



three the lip was too small to measure with the optical microscope mechanism. These 
measurements indicate that as dwell time increases the radial extent of the damage area 
increases, and more material is deposited around the edge of the hole. When wind is 
added to the test, we may observe that debris in the hole is removed during irradiation. 
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V. CONCLUSIONS 



In this thesis, it was established that a new self-defense weapon is needed to 
counter current and future threats to U. S. Navy ships. The FEL was proposed as a 
possible alternative to the inadequate CIWS because of its speed-of-light “bullets”, 
unlimited magazine, tunability, and power. The basic physics of the FEL was described 
through the derivations of the pendulum equation and the optical wave equation. The 
Navy’s Directed Energy Office has already selected the FEL for developmental funding, 
so experimentation is required to determine the direction development should take. This 
thesis describes the first measurements of laser damage from the newly developed 
TJNAF FEL and the results could provide the basis for new directions for directed energy 
weapons design. 

A. SCALING 

The TJNAF FEL, which is capable of several hundred Watts of continuous 
average power, was used to simulate the damage from a MW-class weapon by focusing 
the beam to a smaller spot size. The eventual goal is to develop scaling rules that will 
reliably predict the damage of a large laser without having to bare the enormous cost of 
building the large laser first. The experimental data shows that the scaling concept with 
thermal diffusion calculations is promising. More detailed experiments varying 
wavelength, power, and spot size may be able to produce scaling laws, which would be 
invaluable, for future weapons designers. 
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B. FEL PULSE FORMAT 



The extremely short sub-picosecond pulse length of the FEL beam is a result of 
the electron bunches described in the section on FEL background physics. The TJNAF 
FEL has a unique pulse format with a rapid sequence of short, powerful pulses. The peak 
power in each pulse is about 50 MW lasting for only about one-half picosecond coming 
at a rate of 37 MHz. Other studies have shown that such short pulses may give as much 
as a factor of ten advantage in reduced fluence required to produce damage [Ref. 15]. 

The experiments conducted for this thesis began to collect data to show whether this 
advantage exists, but further experimentation will be required. 

C. FUTURE EXPERIMENTS 

The TJNAF FEL is scheduled for an upgrade to 20 kW of power that will allow 
more flexibility in scaling experiments and further tests of scaling itself. Additional plans 
are for experiments, which include wind passing over the samples, weighing of the 
samples before, and after each run, new wavelengths, changing wavelength during 
irradiation, new pulse formats, and other sample materials. As experimental procedures 
are refined and the amount of data increases, more thorough analysis of the FEL beam 
and comparison to other lasers will become possible. 
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